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Abstract 

We study the production of an isolated large-pT photon as well as the 
production of an isolated prompt photon plus a jet in e“*"e“ collisions. Our 
results are obtained by a NLO Monte Carlo program of partonic event gen¬ 
erator type. We discuss the possibilities to constrain the parton densities in 
the real photon and compare to preliminary OPAL data. 


1 Introduction 


Hard reactions involving real photons, such as the photoproduction of jets K 

BBiB or of large-pr photons [^ |^, are often presented as good tests 

of QCD. Indeed the real photon having a pointlike coupling to a quark of the hard 
subprocess is a good probe of the short distance dynamics. This would particularly 
be the case if one could forget the hadronic component of real photons. In the 
initial state, the incoming real photon can fluctuate into states made of quarks and 
gluons which take part in the hard scattering process. This mechanism, referred to 
as the resolved process, leads to the introduction of quark and gluon distributions 
in real photons. Similarly a hnal state large-pT photon can also be produced by the 
fragmentation of a large-pr parton emerging from the hard scattering process; this 
process is described by introducing the fragmentation function of partons into real 
photons. Therefore the reactions involving real photons are far from being simple 
processes and contain the whole complexity of the pure hadronic collisions. This is 
particularly true for the reaction 77 ^ 7 + X in which the hnal photon is produced 
at large transverse momentum. Each of the photons can have a direct or a resolved 
interaction such that eight different processes contribute to the cross section. 

However this complexity has a positive aspect if the various contributions to the 
cross section can be disentangled by means, for instance, of kinematical cuts. In 
an ideal situation it should be possible to separately observe direct contributions 
which provide unambiguous tests of QCD, resolved contributions leading to the 
measurement of the quark and gluon distributions in the photon and fragmentation 
contributions giving access to the quark and gluon fragmentation functions. These 
fragmentation contributions can be strongly suppressed by using isolation criteria 
HQ which eliminate events with a too large hadronic transverse energy in a cone 
surrounding the photon. For instance the criterion used by the OPAL experiment 
H. that we analyze in this paper, completely suppresses the fragmentation contri- 
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bution and the remaining contributions to the cross section are due to direct and 
resolved incoming photons, the hnal photon being always a direct one. 

One way to disentangle the resolved contribution from the direct one rests on 
the measurement of the momentum of the parton entering the hard subprocess [Q 
(a measurement which requires the observation of the whole hnal state). When 
the ratio of this momentum to the incoming photon momentum is close to one, one 
expects the direct contribution to be dominant (the photon momentum is not shared 
among the quarks and gluons); for a ratio smaller than one the resolved contribution 
is the larger one. One interesting aspect of the 7 + 7 —> 7 + X reaction at LEP II 
energies is that the dominant contributions to the cross section contain at least one 
resolved photon in the initial state. The direct-direct correction is negligible, and 
the observation of the 7 -|- 7 —>• 7 -|- X cross section opens the possibility to measure 
the parton distributions in the photon. 

Another interesting point of the reaction 7-1-7^7-l-Xat LEP II energies is the 
small contribution of the NLO corrections. As we shall see below, these corrections 
represent - for the scales set equal to pt^ - about 10 % of the cross section. This 
result shows that the QCD corrections are well under control. Moreover the scale 
dependence is very weak and therefore the theoretical predictions are unambiguous. 

The possibility to constrain the parton distributions in the photon, especially the 
rather poorly known gluon distribution, via the process 7 -f 7 —>■ 7 -|- X has already 
been pointed out in |jT^, While |T^ is a leading order calculation, in the 
calculation is done at next-to-leading order, but only the fully inclusive cross section 
is calculated, without the possibility to consider isolated photons or a photon -|- jet 
final state. 

The plan of the paper is the following. Section 2 is devoted to the theoreti¬ 
cal framework and to the discussion of predictions for the inclusive cross sections 
da/dpT^ and da/dp^. In Section 3 we study the reaction 7 -|- 7 —> 7 -|- jet + X, 
emphasizing the difference to the inclusive case. Comparisons with preliminary re- 
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suits of the OPAL experiment are performed in Section 4 . Section 5 contains the 
conclusion. 


2 Theoretical Framework and Predictions 

2.1 General Setting 

As the general framework of the calculation has already been described in detail 
,|,i, we will give only a brief overview on the method here. 

In reactions, the electrons can act like a source of quasi-real photons whose 
spectrum can be described by the Weizsacker-Williams formula 
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This approximation is valid as long as Qmax/PT 1 where pt is a large scale 
characterizing the hard subprocess, for instance the transverse momentum of the 
hnal photon. In this section we use this approximation, postponing to Section 4 a 
discussion of its validity when comparing our predictions with OPAL data. 

The quasi-real photon then either takes part directly in the hard scattering pro¬ 
cess, or it acts as a composite object, being a source of partons which take part in 
the hard subprocess. The latter mechanism is referred to as resolved process and is 
parametrized by the photon structure functions Fa/^{x^, Q"^). Thus the distribution 
of partons in the electron is given by the convolution 


Faieixe, M) = f dy dx^ f^{y) Fa/-y{x^, M) 5 {x^y - Xe) ( 2 ) 

^ 0 

where in the “direct” case Fa/^{x^, M) = (5a7(5(l — x^). 

Similarly, a high-pT photon in the final state can either originate directly from 
the hard scattering process or it can be produced by the fragmentation of a hard 
parton emerging from the hard scattering process. However, as discussed in the 
introduction, this fragmentation contribution can strongly be suppressed by using 


3 








an isolation criterion. It tnrns ont that this isolation is also nsefnl at the experimen¬ 
tal level in order to single ont the prompt photon events from the backgronnd of 
secondary photons prodnced by the decays of light mesons. The isolation criterion 
nsed by the OPAL collaboration is the one proposed by Frixione which leads to 
a complete snppression of the fragmentation contribntions (contribntions involving 
the fragmentation fnnction of parton a into a photon). At the parton 

level we have the following constraints : for each parton i, the distance 


Rh = V i<Pi - (t>iY + iVi - Vi? ( 3 ) 

to the photon is compnted in (0, rj) space, where 0 and r] are the azimnthal angle 
and psendorapidity, respectively. A photon is kept if the condition 

^ - Ri,^) < 0.2 ■ Et,^ I ~ , forain<i? ( 4 ) 

i cos^irj 

is fnlhlled, where Ej-^i is the transverse energy of the parton, 0 is the step 
fnnction, which ensnres that only particles in the cone with apertnre 6 contribnte 
to the snm, and the cone radins is i? = 1. 

As a resnlt, we only have to consider contribntions with a direct hnal photon. 
These can be classihed by the three prodnction mechanisms illustrated in Fig. 1 : 1 ) 
direct-direct, 2 ) single resolved, 3 ) double resolved. Note that by ’’single resolved” 
we denote the sum of direct-resolved and resolved-direct contributions. 

We implemented the full set of next-to-leading order corrections to the single resolved 
and double resolved processes. Note that the direct-direct process is by itself a higher 
order correction; there is no Born contribution involving a 2 —> 2 subprocess in this 
reaction. The isolation criterion does not modify the leading order contributions for 
the direct-resolved and resolved-resolved case shown in Fig. |l|. But at NLO a parton 
can be close to the photon in the (0, rj) plane and the cross section is modihed by 
the constraint (|^). 
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Figure 1: Examples of contributing processes. 


In the calculation of the higher order corrections a phase-space slicing method 
is used which allows to isolate and analytically calculate the divergent infrared and 
collinear contributions. The non-divergent parts of the 2 —>■ 3 contributions are 
calculated numerically. A partonic event generator (allowing negative weights) has 
been built on these bases. The interested reader may hnd a detailed description of 


this generator in ref. [16 


2.2 Predictions 

Let us now study the NLO cross section of the —*• e^e~^X reaction obtained 

in the above theoretical framework and by using the OPAL kinematical conditions. 
Thus we use \/Ag+e- = 196.6 GeV, — 10 GeV^ and 0 < ?/ < 1. We also use 
the AFG parton distributions and the value A^^ = 300 MeV. (We work with 4 
flavours). The factorization scale M and the renormalization scale /i are set equal 
to pt^. The NLO cross sections da/dpr^ associated with the different production 
mechanisms, integrated in the rapidity range —1 < 1, are displayed in Fig. 2. 

For the scales p = M = px^ used here the single resolved contribution is larger 
than the double resolved contribution by a factor 3 to 6 which increases when pt^ 
increases. The single resolved contribution constitutes about 70 % of the total cross 
















Figure 2: The NLO single resolved, double resolved and direct cross sections 
/dpT^. 

section. 

The ratio of the HO corrections to the NLO (= LO + HO) cross sections is 
displayed in Fig. |[ The HO corrections are small for the single resolved cross 
section (HO/(LO + HO) ~ 0.1)) and large for the resolved-resolved case (HO/(LO 
-|- HO) ~ 0.5). The direct-direct HO contribution is smaller than the single resolved 
HO contribution; it becomes negative at large pt^ {pt^ > 10 GeV). On this occasion 
we remind the reader that the individual contributions to the cross section are 
not physical. They depend on the factorization and renormalization scheme (here 
the MS scheme). For instance the direct-direct HO contribution comes from the 
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Figure 3: Ratios HO/(LO + HO) for the e’''e —*• e'''e 7 X cross section, calculated 
with the scale /r = M = as a function of px^. 

7 + 7—^7 + g + g cross section from which the initial state collinear divergences 
have been subtracted in the MF-scheme. The hnite remainders of the subtracted 
singularities depend strongly on the factorization scale M, this dependence being 
compensated only when adding the single resolved contributions. 

Let us turn to the question of the stability of the predictions with respect to 
the choice of the renormalization and factorization scales. In Fig. ^ we display the 
variation of the cross section da/dpx^ calculated at px^ = 6 GeV as a function of 
the parameter C which dehnes the scale choice, p, = M = C px^ resp. M = C px^ ; 
the cross sections are normalized to one at G = 1. The increase of the leading 
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Figure 4: The dependence of /dpT^ at px^ = 6 GeV on the renormalization 

and factorization scales as a function of C. Solid line: p = M = Cpx (NLO), 
dotted line: p = M = Cpx^ (LO), dashed line: p hxed to px^, M varied between 
PxJ^ and 2,px^ (NLO). 


ordei|3 cross section with increasing scales reflects the fact that - in contrast to the 
hadronic structure functions which show the well known scaling violations ~ the 
photon structure functions grow uniformly like log(Q^/A^) for large since the 
pointlike piece dominates at large Q^. 

We note that the variations of the NLO cross section are smaller than 5 %, when 
the values of the scales M and p are changed simultaneously by a factor 9. At the 

same time, each individual cross section strongly varies. For instance the direct- 

*What we denote by ’’leading order” only refers to the matrix elements. For the parton distribu¬ 
tions we always use the NLL fits, and for as we use an exact solution of the two-loop renormalization 
group equation, and not an expansion in log/r/A. 
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direct contribution varies from 7.85 xlO“^pb to —6.67 x 10“^ pb whereas the total 
cross section is equal to 23.82 xl0“^pb at C = 1. If we vary the factorization scale 
M only, keeping /i hxed to , the variations of the NLO cross section are of course 
larger since cancellations of logp'^^/M'^ and log/r^/pf.^ terms are spoiled, but they 
are still only of the order of 20% for a variation of by almost two orders of 
magnitude, as can be seen from Fig. 

To conclude this short study of the cross section da/dpT^ we can say that the NLO 
calculation is well under control; the HO corrections are small and the sensitivity of 
the cross section to the scales is negligible. Therefore the reaction —>■ e'^e~'yX 

offers interesting tests of QCD and should allow to put constraints on the parton 
distributions in the real photon. 

The cross section da/dp-y is a good observable to determine Fa/^{x, M), because 
the value of x is closely related to the value of As a consequence the shape 
of da/dp.y gives indications on the x-dependence of the parton distributions in the 
photon. The cross sections da/dp-y are shown in Fig. |]for px^ > 3 GeV. Obviously 
the result is symmetric with respect to = 0 within the statistical errors. The box 
contribution pg pg to the single resolved part is of the order of 5% of the total 
single resolved contribution, and of the same order of magnitude as the direct-direct 
contribution. However, in view of other uncertainties of the order of 5 % (see Section 
4), it has not been included in the results shown. The box correction gg —>■ pg to 
the double resolved part is one order of magnitude smaller than the 'pg —>■ pg box 
and therefore completely negligible. 

Large values of p^ enhance the direct-resolved contribution (where the incoming 
direct photon goes towards positive rapidity). The resolved-resolved contribution is 
quite flat and the direct-direct one almost negligible (at the scales p = M = px^)- 
A value of larger than 3 GeV would enhance the single resolved cross section 
with respect to the double resolved one. 
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Figure 5: The cross section /dr]^ as a function of r]^ for pt^ > 3 GeV. 





















3 Photon—jet cross section 


In the introduction, we discussed the possibility to enhance a resolved contri¬ 
bution with respect to a direct one by hxing the momentum of the partons entering 
the hard subprocess. Following the OPAL collaboration, we dehne 


^ _ {Ejet ± + {E^±p2) _ + p?’ 

“ E {E±p,) ~ ^ ^ 

hady'y 

where are the fractional momenta of the quasi-real initial photons oriented to¬ 
wards the positive and negative z-axis, y = E^/E^. The conventions are such that 
particle one travels towards the positive z-direction, such that e.g. direct-resolved 
means that the direct photon - with momentum fraction - moves to positive 
z whereas the resolved photon (with momentum fraction x~) moves to negative . 2 . 
The variables x^ dehned by eq. are exactly the longitudinal momentum fractions 
(with respect to the photon momenta) of the partons entering the hard subprocess 
if the latter is a 2 —2 process. In the case of a 2 3 subprocess, these momentum 

fractions are no more hxed by the observation of the photon and a jet since there 
are conhgurations with a third (unobserved) particle in the hnal state such that the 
true value of x^ in Ea/^{x^, M) (see eq. (@)) can be larger than Nevertheless 
remain useful variables to constrain the parton momenta. 

The dehnition (^ implies to measure the jet rapidity and transverse energy. In 
ref. we dehned another variable which does not require the measurement of p^*. 


^LL = 


y^ \/ Se-e+ 


( 6 ) 


This is certainly an advantage with respect to expression (^, because the mea¬ 
surement of a jet transverse energy at small p^^* (p^ * ~ 5 GeV) can be very difficult 
and inaccurate ; it is delicate to disentangle the jet energy from the underlying event 
energy. Moreover the variables (^ lead to smoother theoretical distributions in the 
regions ~ 1, which may be useful if one wants to compare theory with data for 
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these values of These points have been discussed in more detail in ref. to 

which we refer the interested reader. 

Before using these variables, let us consider the behaviour of the photon-jet cross 
section with respect to the photon and jet rapidities. In order to enhance the direct- 
resolved contribution, we choose a large value of rjjet, 1 < rjjet < 2, and we vary 
rjy. The jets are dehned in agreement with the cone algorithm used by the OPAL 
collaboration []^ , with a cone radius R = 1 and a cut on the minimum value of 
that we £x equal to 3 GeV |^. We also £x = 5 GeV. Our results are shown in 
Fig. I and they can be compared to those displayed in Fig. ||. 
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Figure 6: The cross sections /di]^ at pt^ = 5 GeV and 1 < rjjet < 2 for two 

different scale choices. 
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We recall that, as always, the individual contributions to the total cross section 
have no physical meaning. Their relative importances are scale dependent. However, 
choosing e.g. the scales /i = M = does not modify the fact that the direct- 

resolved contribution dominates in the range ri^>— 0.7. Fig. ^ shows the enhanced 
direct-resolved contribution, corresponding to small values of and large values 
of 2 :^^) in expressions (|^) or (||), for the two scale choices p = M = pt^ and 
p = M = Pt^/2. Hence, by working at large photon and jet rapidities, we can 
enhance the direct-resolved contribution and study its behaviour as a function of 
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Figure 7: /dx^ with px^ = 5GeV, 1 < pjet < 2 and 0.5 < <1. 

In Fig. 0 we display the distribution da/dx^^ obtained with px^ = 5 GeV, 1 < 
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Vjet < 2 and 0.5 < ^77 < 1. The contributions from direct-direct and resolved- 
direct have been summed up in Fig. ^ since at ~ 1 there are large cancellations 
between these two contributions due to the finite remainders of the initial state 


collinear singularities, as has been explained in Section |2.1| . This is illustrated in 
Fig. ^(a). Of course there are also cancellations between the direct-resolved and the 
direct-direct part, but those are of major importance in the x'^ spectrum for x'^ ~ 1 , 
affecting the x~ distribution only marginally. 

Thus we can conclude from Fig. ^ that in the domain xJ^j^ ^0.8 the cross section 
/dxli is quite sensitive to the parton distributions Fa/^{x~ ~ x2l,Pt^) in 
the photon. 

Until now we have used the variable xll to study the various contributions to 
the cross section da/dx. The use of x^ as defined in eq. (|^) leads to a similar result, 
except around x ~ 1 where the cross section da/dx^i has a smoother behaviour, as 
illustrated in Fig. ||(b). 

Fig. H shows that the contribution from the subprocess 7 -I- 5 ' —> q + + q 

contributing to the direct-resolved part of da/dr]^ is very small, and even negative 
for ?77 > 0 , such that the full direct-resolved contribution is slightly lower than the 
one where g'^ has been set to zero. Thus the sensitivity to the gluon distribution in 
the photon as shown in Fig. ^(a) is basically due to the double resolved part in the 
present study. 
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(a) cancellations between resolved-direct and direct-direct subprocesses 

(b) comparison between and /dx~ (for the sum of all subpro¬ 

cesses). The distribution in xll shows a much smoother behaviour at a; ~ 1. 
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4 Comparison with OPAL data 


Before comparing theoretical predictions and data we have to assess the validity 
of the Weizsacker-Williams (WW) approximation embedded in formula (|^). This 
approximation is obviously correct as long as the photon virtuality is negligible with 
respect to the scale involved in the hard subprocess, namely pt^. But in OPAL 
Qmax — 10 GeV^ is of the order of the subprocess scale, pt^ > 3 GeV, and higher 
twist contributions proportional to Q^ax/Pr^ ^ priori negligible. 

However we have to note that the Weizsacker-Williams formula is dominated by 
the term ~ /ml) which is very large (~ 27) because of the small mass 

of the electron; therefore the WW approximation should not be too bad, even if 
Q'LaxIPT^ - 1 - 

Moreover the large virtuality of the incoming photon questions the use of parton 
distributions in real photons. To estimate the importance of the corrections due to 
the photon virtuality, let us consider the perturbative quark distribution in a real 
photon in the Born approximation 


q,{x,M) = 3-{x +(l-^))4 ^ • (7) 

The lower limit is the boundary below which a perturbative approach is no 

more valid. The value Ql ~ m^ ~ 0.5 GeV^ leads to a photon structure function 

F^{ Q"^) = X [Q'yi.x, Q)+Q-y{x, Q)] in agreement with DIS 7*7 experiments 

g 

For a virtual photon, Ql must be replaced by max{Q‘^, Ql) and the convolution with 
the spectrum of the photon emitted by the electron takes the form 


[171 


I = 


,2 

max 


dQ^ 


rM^ 


dk‘^ 


( 8 ) 


J max{Q'^ ,Q'^) kP' 

in which we drop all x-dependence. If > Ql, expression (||) can be written as 
the sum of the real photon expression and a correction 
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I = 


2 

max 


dQ^ dk^ 1 2 QL. 

Q2 Jq 2 k^ 2 Ql 


(9) 


For Q'^ax = 10 GeV^, Ql = 0.5 GeV^ and = 25 GeV^, the correction is 

of the order of 5 %. A similar estimation can be done for the non-perturbative part 
of the photon strnctnre fnnction with a correction of the order of 15 % which is also 
negative . 

The conclusion of this study of the effects of the initial photon virtualities is that 
the WW approximation is valid in the OPAL kinematical range up to corrections of 
a few tens of percents. 

Now we compare our results to OPAL preliminary data [jl^. It has to be noted 
that in the experimental analysis the aim was to keep only events from single and 
double resolved processes. Therefore a component identified as stemming mainly 
from direct-direct events (denoted by ”£nal state radiation” (FSR) in [^) has 
been subtracted from the data. In a NLO calculation the final state radiation 
associated with the subprocess 77 —> 'yqq is a genuine contribution to the higher 
order corrections and thus cannot be removed from the total cross section. It can 


only be suppressed by an isolation criterion as discussed in Section PT| . In the present 
study the isolated direct-direct contribution is small, such that it should not disturb 
the comparison with the preliminary OPAL data which have large experimental 
errors. 

In Fig. 1^ the cross section da'^^^/dpT^ is shown. One observes that NLO theory 
is somewhat higher than the data, in particular at low px ^, but a correction term due 
to the rather large value of Q^ax would decrease the theory prediction (see eq. (|])), 
and including the FSR in the experimental analysis would increase the data. On 
the other hand, it is hard to estimate if higher twist effects (which should be non- 
negligible only at low pxQ would increase or decrease the theory prediction in the 
first pT-y-hin. 

In short, taking into account a theoretical uncertainty of about 10%, the agreement 
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Figure 10: Comparison to preliminary OPAL data for /dpr^ integrated over 

photon rapidities in the range —l<ri^< 1, at the scales p = M = pt^. 


with the data is completely satisfactory. The comparison to the rapidity distribution 
da'^+^/d\py\ is displayed in Fig. 0. 

In [0] it is shown that PYTHIA reproduces the shapes of the distributions well, 
but has to be scaled up by a factor 1.85 to be consistent with the data. One reason 
for this difference in normalisation might be the use of different parton distribution 
functions for the photon (to obtain the PYTHIA result, SAS-ID - which are 
LO hts - have been used). 

We further note that the data are also compatible with our ’’leading order” result. 
This does not come as a surprise since the higher order corrections are small and 
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Figure 11: Comparison to preliminary OPAL data for /d\iri.y\, at the scales 

li = M = pt^. 

the preliminary data fluctuate a lot, but we recall that by ’’leading order” we only 
refer to the matrix elements and not to the parton distribution functions and o;*. 
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5 Conclusions 


We have presented a complete NLO program for the calculation of the reactions 
e’''e“ —»■ e’''e“ 7 X and —*• e’''e“ 7 + jet + X, where the prompt photon is 

isolated using the criterion of S. Frixione []^. The higher order corrections turned 
out to be of the order of 20% and the scale dependence is very weak: The total cross 
section varies by less than 5% if the scales are varied between p, = M = Pt^/3 and 
^Pt^- 

We discussed the possibility to constrain the parton distributions in the photon. 
For the photon-jet cross section, we have shown that one can enhance the direct- 
resolved component of the cross section by restricting the photon and jet rapidities 
to large values. 

Further, we discussed the validity of the Weizsacker-Williams approximation 
for a comparison to OPAL data with a rather large maximal photon virtuality of 
Qmax = lOGeV^. The theoretical uncertainties introduced by using the parton 
distributions of real photons are also estimated in this context. 

We compared the photon px and rapidity distribution for the 77 —>■ 7 -|- X cross 
section to preliminary OPAL data |]^. In view of the large experimental errors, the 
agreement is good. 
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